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Abstract. 

Magnet ometry and magnetic imaging with nitrogen- vacancy (NV) defects in diamond 
rely on the optical detection of electron spin resonance (ESR). However, this technique 
is inherently limited to magnetic fields that are weak enough to avoid electron spin 
mixing. Here we focus on the high off-axis magnetic field regime for which spin 
mixing alters the NV defect spin dynamics. We first study in a quantitative manner 
the dependence of the NV defect optical properties on the magnetic field vector 
B. Magnetic-field-dependent time-resolved photoluminescence (PL) measurements are 
compared to a seven-level model of the NV defect that accounts for field-induced spin 
mixing. The model reproduces the decreases in (i) ESR contrast, (ii) PL intensity 
and (iii) excited level lifetime with an increasing off-axis magnetic field. We next 
demonstrate that those effects can be used to perform all-optical imaging of the 
magnetic field component \B±\ orthogonal to the NV defect axis. Using a scanning 
NV defect microscope, we map the stray field of a magnetic hard disk through both 
PL and fluorescence lifetime imaging. This all-optical method for high magnetic field 
imaging at the nanoscale might be of interest in the field of nanomagnetism, where 
samples producing fields in excess of several tens of milliteslas are typically found. 
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1. Introduction 

The nitrogen-vacancy (NV) defect in diamond is a solid-state quantum system that 
has been extensively studied over the last years |1|. Potential applications include 
quantum information processing (2H5], imaging in life science |6|, hybrid quantum 



systems [7j-[9j and magnetic sensing and imaging 10-13 . For the latter application, the 
magnetic field is evaluated by measuring the Zeeman shifts of the NV defect electron 
spin sublevels through the optical detection of electron spin resonance (ESR). The 
main advantage of NV-based magnetometry is the possible combination of atomic- 
scale spatial resolution with high magnetic field sensitivity - below 10 nT/Hz 1 / 2 [l4] - 
even under ambient conditions. Magnetic imaging with diffraction-limited micrometer 
resolution has been demonstrated using an ensemble of NV defects [15-17 while recent 



experiments have reached the nanoscale with a single NV defect coupled to a scanning 
probe microscope [l3[[l8|[l9) . 

However, ESR-based magnetometry with NV defects is intrinsically limited to 
magnetic fields with an amplitude and an orientation such that the electron spin 
quantization axis remains fixed by the NV defect axis itself. Indeed, any significant spin 
mixing induced by an off-axis magnetic field rapidly reduces the contrast of optically- 
detected ESR spectra because optically induced spin polarization and spin dependent 
photoluminescence (PL) of the NV defect become inefficient. Besides a decreased ESR 
contrast, the PL intensity as well as the effective excited level lifetime are observed to 
decrease with an increasing off-axis magnetic field (22j[23). It was recently shown that 
this property can be used as a resource to perform all-optical magnetic field mapping 



with a scanning NV defect [18, 19 . The purpose of the present work is to study in a 
quantitative manner the dependence of the NV defect optical properties on the magnetic 
field vector B, and to discuss the implications for magnetic field imaging. The paper is 
organized as follows. In Section [2j we first investigate the dynamics of optically-pumped 
NV defects as a function of B, focusing on (i) the ESR contrast, (ii) the PL intensity 
and (iii) the excited level lifetime. A simple seven-level model of the NV defect is 
developed and compared to magnetic-field-dependent time-resolved PL measurements. 
This method provides a simple way to extract all relevant photophysical parameters of 
the NV defect and compute its optical response as a function of the B field. The results 
of the calculation are then compared to single site measurements and discussed in the 
context of magnetic field imaging. In Section [3j we report on experiments of all-optical 
magnetic imaging in the high off-axis magnetic field regime using a scanning NV defect 
microscope. In particular, fluorescence lifetime imaging (FLIM) capabilities allows us 
to record not only PL images but also lifetime ones, a feature that could reveal useful 
beyond magnetic imaging, e.g. for mapping the local density of electromagnetic states 



(LDOS) of photonic nanostructures [20,21 . 
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2. Photodynamics of single NV defects in a static magnetic field 

2.1. Model 

The negatively charged NV defect in diamond consists of a substitutional nitrogen atom 
(N) associated with a vacancy (V) in an adjacent lattice site of the diamond matrix, 
giving a defect with symmetry (Fig. [j](a)). Its ground level is a spin triplet 3 A 2) 
whose degeneracy is lifted by spin-spin interaction into a singlet state of spin projection 
m s = and a doublet m s — ±1, separated by D gs = 2.87 GHz in the absence of magnetic 
field (24j. Here, m s denotes the spin projection along the NV defect axis z (Fig. [jja)). 
When a magnetic field B is applied to the defect, the ground level spin Hamiltonian 
reads 

n gs = hD gs s 2 z + 9f i B B -s = n z gs + (1) 

with % z gs = hD gs S 2 + gjiBB z S z and 7i gs = gjJLB{B x S x + B y S y ), where h is the Planck 
constant, \ib is the Bohr magneton, and g « 2 is the electron ^-factor. We note 
that strain-induced splitting and hyperfine interaction with nearby nuclear spins are 
neglected in Eq. ([!]). In zero field, the quantization axis is thus fixed by the NV defect 
axis and the eigenstates are labelled |2°) and |3°), corresponding to spin projections 
m s = 0, m s = —1 and m s +1, respectively (Fig. [jjb)). 

The defect can be optically excited through a dipole- allowed transition to a 
excited level, which is also a spin triplet. Besides, the 3 E excited level is an orbital 



doublet which is averaged at room temperature 25 , 26 leading to a zero-field splitting 
D es = 1.42 GHz with the same quantization axis and similar gyromagnetic ratio as in 
the ground level [27j[28]. The excited level spin Hamiltonian H es is thus simply given 
by Eq. while replacing D gs by D es . The excited level eigenstates in zero field are 
labelled |4°), |5°) and |6°), corresponding to spin projections m s = 0, m s = — 1 and 
m s = +1, respectively (Fig. [j](b)). Once optically excited in the 3 E level, the NV 
defect can relax either through the same radiative transition producing a broadband 
red photoluminescence (PL), or through a secondary path involving non radiative 
intersystem crossing (ISC) to singlet states. Recent experiments have identified two 



singlet states i^E^Ai) [29, 30 , whereas theoretical studies predict the presence of a 



third singlet state between the ground and excited triplet levels 31 . In the present 
work, the singlet states are summarized into a single 'metastable' level labelled |7). 
Spin-dependent ISC from the 3 E level to the 'metastable' level is responsible for efficient 
spin-polarization in the m s = spin sublevel through optical pumping as well as for 
spin-dependent PL of the NV defect. These two properties enable the detection of 
electron spin resonance (ESR) on a single NV defect by optical means [32] . 

Within this seven-level model of the NV defect, the zero-field transition rate from 
\i°) to is denoted by k^. We consider that optical transitions are purely spin- 
conserving, and that the radiative relaxation rate is spin-independent, i.e. k\ x = k® 2 = 
fcg 3 = k®. This is reasonable since Robledo et al. recently determined that spin-flip 
radiative transition rates are at most 2% of their spin-conserving counterpart [33] . On 
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Figure 1. (a) Atomic structure of the NV defect in diamond. A magnetic field B 
is applied with an angle with respect to the NV defect axis z. (b,c) Seven-level 
energy level structure considered for the NV defect at room temperature. The ground 
( 3 A2) and excited ( 3 E) spin triplets are linked by spin conserving radiative transitions 
(solid lines) while the singlet states, merged into state |7), provide a non-radiative 
relaxation path (dashed lines). All notations are defined in the main text, (b) At 
B = 0, five intrinsic parameters {fcj, &4 7 , fc{- 7 , fc 71 , k 72 } are used to describe the NV 
defect photodynamics. (c) For a given applied magnetic field £?, all transition rates kij 
are likely to be non vanishing owing to electron spin mixing both in the ground level 
and in the excited level. 



the other hand, the optical pumping rates from the ground to the excited level are 
proportional to the corresponding relaxation rates through k^ = /3k^ for j = 1,2,3 
and i — 4, 5, 6, where f3 is the optical pumping parameter. Lastly, ISC transition rates 
at zero- field only depend on the absolute value of m S) i.e. k® 7 = k® 7 and k 72 = ^73- 
Following these assumptions, there are five intrinsic parameters left in the system, 
namely k® 7) k® 7) k 71 and k 72) and an extrinsic parameter f3 linked to the optical 
pumping power (Fig. [jjb)). 

When a static magnetic field B is applied to the NV defect (Fig. [jjc)), the seven 
eigenstates {\i)} of the system can be expressed as linear combinations of the zero-field 
eigenstates 

K) = ^^-(B)b°), (2) 

where the coefficients {a^-(B)} are numerically computed by using the expressions of 
Hgs and Hes, and given that \7) = |7°). The new transition rates {fcy(B)} are then 
related to the zero-field transition rates through the transformation 

7 7 

M B ) = ^J2\ a ip\ 2 \ a 3q\ 2k pq ' ( 3 ) 

p=l q =l 
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which represents a statistical averaging using the weights of the initial and the final 
state of the transition. The lifetime of each eigenstate in the excited level is then given 
by n = l/(X)J=i hj) with i = 4, 5, 6. 

With the aim of studying the time-resolved optical response of a single NV defect 
placed in a static magnetic field, we consider a pulsed optical excitation with a pulse 
duration St and a repetition period T such that T » Tj » St. After optical excitation 
(t — 0), each excited state population rii decays exponentially as r^(t) = r^(0)e -t / Ti 
once averaged over many excitation cycles. The rate of emitted photons 7t^(i,B) from 
state \i) is the radiative part of the total decay and reads 7t^(t, B) = J^ =1 %r^(i). The 
detected time-resolved PL signal 1Z(t, B) is thus finally obtained by summing over the 
three excited states 

6 3 

K(t,B) = %"i(0)e-* /T< , (4) 

i=4 j=l 

where rj is the collection efficiency. In this equation, the magnetic field dependent 
transition rates {%} are given by Eq. ^ and the set of coefficients {^(0)} are inferred 
by solving the classical rate equations of the system 

— -- = ^^(kjiTij — kijfii) , (5) 

while considering a closed seven-level model, i.e. Yli=i n i = 1? an d the periodic condition 
rii(t) = rii(t + X). The detected time-resolved PL signal 1Z(t, B) is therefore given by 
a tri-exponential decay, where the amplitude and the lifetime of each component is 
expressed as a function of the zero-field transition rates {k®, k® 7 , k® 7 , k 71 , k 72 } and the 
optical pumping parameter f3. We expect that 77,4 3> (77,5, uq) at low magnetic field owing 
to efficient spin polarization of the NV defect in m s = by optical pumping, while those 
populations are balanced in the limit of strong transverse magnetic fields. 

Finally, we infer the contrast C of optically-detected ESR spectra for a single NV 
defect placed in a static magnetic field. For that purpose, non-vanishing transition 
rates fc 12 = k 2 i (resp. fc 13 = fc 31 ) are added in the model to account for a microwave 
field in resonance with the |1) <R> |2) (resp. |1) |3)) electron spin transition in the 
ground level. In the following ESR spectra are recorded with continuous-wave (CW) 
optical excitation. In this case, the mean PL rate is given by 7£(B) = r\ Y^i=A Y^j=i ^ikij 
where the averaged populations {n^} are inferred by solving the classical rate equations 
(Eq. ([5])) at the steady state. In this framework, the ESR contrast of the |1) |2) 
transition is defined by 

C(B) = — Wh^) — ' (6) 

where TZ(ku = 0) (resp. 1Z(ki 2 )) denotes the NV defect PL rate without applying the 
microwave field (resp. with a resonant microwave field). A similar expression holds for 
the |1) |3) transition. 
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2.2. Experimental setup 

Individual NV defects hosted in a high-purity type Ha diamond crystal (Element Six) are 
optically isolated at room temperature using a scanning confocal microscope. Optical 
excitation at 532 nm wavelength is provided either by a CW laser or by a pulsed laser 
with a pulse duration St « 60 ps and a \jT — 10 MHz repetition rate (PicoQuant, LDH- 
P-FA-530B). The detection system comprises a confocal arrangement, an avalanche 
photodiode working in the single-photon counting regime (Perkin-Elmer, SPCM-AQR- 
14), and a time-correlated single-photon counting module with a bin size of 512 ps 
(PicoQuant, PicoHarp 300). An electromagnet is used to apply a static magnetic field 
with controlled amplitude, while ESR transitions in the ground level are driven with a 
microwave field applied through a copper microwire directly spanned on the diamond 
surface. 



2.3. Results and discussion 

For each studied NV defect, the angle 9 between its symmetry axis (z) and the magnetic 
field B was first measured by using ESR spectroscopy (see Tab. For that purpose, 
optically-detected ESR spectra were recorded by sweeping the frequency of a microwave 
field while monitoring the NV defect PL intensity. Owing to spin dependent PL, two dips 
can be observed in the PL signal, corresponding to electron spin transitions |1) <H> |2) and 
|1) 1 3) (Fig. [2](a)). Measuring the electron spin transition frequencies as a function 
of the magnetic field amplitude B = ||B|| allows us to extract the angle 9 through data 
fitting with the eigenenergies of the ground level spin Hamiltonian H gs (Fig. gb)) [l3|. 

The ESR spectra depicted in Figure [2](a) can also be used to estimate the ESR 
contrast C as a function of the magnetic field amplitude (Fig. [2](c)). If the magnetic 
field is such that H| ses 3> 7^ >es (see Eq. ([!)), the quantization axis is fixed by the 
NV defect axis, and a high ESR contrast is observed. Conversely, when the condition 
T~L z gs es 3> 1-Lg S es is not fulfilled, the quantization axis is rather determined by the applied 
magnetic field and m s is no longer a good quantum number. The eigenstates of the 
spin Hamiltonian are then given by superpositions of the m s = and m s = ±1 spin 
sublevels, both in the ground and excited levels (see Eq. ([2])). As a result, optically- 
induced spin polarization and spin-dependent PL of the NV defect become inefficient, 
and the contrast of optically detected ESR vanishes, as shown in Figure [2](c). Magnetic 
field imaging through optically-detected ESR is therefore inefficient in the regime of 
'high' off-axis magnetic field. Nevertheless, the NV defect optical response can be used 
to extract information about the magnetic field in this regime, as explained below. 



The photophysical parameters of the NV defect were first measured in order to com- 
pute its optical response as a function of the B field. For that purpose, time-resolved PL 
traces 7£(£, B) were recorded while increasing the magnetic field amplitude B for a given 
angle as shown in Figure [3](a). Fitting the whole set of data with Eq. Q allows us 
to extract the zero-field transition rates {k®, k® 7 , k® 7) k 7l) k 72 } as well as the normaliza- 
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Microwave frequency (GHz) g (mT) 

Figure 2. (a)-Optically detected ESR spectra recorded for a single NV defect (NV D) 
at various amplitudes B of the magnetic field. Solid lines are data fitting with 
Lorentzian functions. (b)-ESR frequencies as a function of B. Solid lines are data 
fitting with the eigenenergies of 7~L gs . The angle deduced from the fit is = 74 ± 1°. 
(c)- ESR contrast as a function of B. Markers are data corresponding to the |1) o \2) 
electron spin transition while the solid line is a fit using Eq. ([6| of the seven-level 
model. Here the only free parameters are the optical pumping parameter j3 and the 
ESR transition rate k\2. All other parameters of the model are obtained through time- 
resolved PL measurements as summarized in Table 1. For these experiments the NV 
defect was excited with a CW laser power of 100 /iW, leading to a detected PL intensity 
of 5 x 10 4 counts. s _1 at zero field. The evolution of the normalized PL intensity of 
NV D as a function of the magnetic field strength is shown in Fig. pVc) . 



tion coefficient r) and optical pumping parameter f3 (Fig. |3](b)). Horizontal and vertical 
linecuts of the matrix 7£(i, B) are plotted in Figure [3](c) and (d), respectively, showing 
good agreement between the experimental data and the model. We emphasize the fact 
that such a tri-exponential fit can be reliably applied because all traces are fitted at 
once. We thus benefit from the correlations between traces at different magnetic field 
amplitudes, whereas each trace taken individually could not be fitted satisfactorily with 
a tri-exponential function. The transition rates extracted from the fitting procedure 
are summarized in Table [T] for four different NV defects oriented with different angles 
9 with respect to the magnetic field. We note that these values are in agreement with 
those reported in Ref. (33), which were obtained using a different approach based on 
multi-pulse excitation techniques combined with coherent electron spin manipulation. 
As expected, we measure a spin-dependent ISC rate k® 7 3> k® 7 , which is responsible 
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Detection time (ns) B (mT) 

Figure 3. (a) Time-resolved PL traces 7£(£, B) measured for NV E as a function of 
the magnetic field amplitude B. PL counts, corresponding to a 15 s integration time, 
are encoded in a logarithmic color scale, (b) Fit to the data as explained in the main 
text. The reduced chi-square error of the fit is Xred = 1-15- The parameters extracted 
from the fit are summarized in Table [I] (c) Linecuts corresponding to horizontal 
dashed lines in (a) and (b) showing several decay traces measured for various values 
of B (markers) together with the fit (solid lines). At zero-field (red line), the relative 
amplitudes of the three exponential decays are 0.72, 0.14 and 0.14, corresponding to 
decay times T4 = 13.7 ns, T5 = 8.6 ns and tq = 8.6 ns, respectively. At the highest 
field (black line) they are found to be 0.34, 0.34 and 0.32, with decay times T4 = 11.0 
ns, T5 = 9.4 ns and tq = 9.1 ns, respectively, (d) Linecuts corresponding to vertical 
dashed lines in (a) and (b) showing measured PL counts as a function of B at various 
detection times (markers) together with the fit (solid lines). 

for optically-induced spin polarization and spin-dependent PL response of the NV de- 
fect. Furthermore, our measurements seem to confirm that the metastable state decays 
roughly as often in the m s = state as in the m s = ±1 state, i.e. k® r ~ k® 2 [33], in 
contradiction with the initial belief in a coupling to m s = alone [24] . 

Having completely characterized the NV defect and extracted its photophysical pa- 
rameters, we now discuss the dependence of the NV defect optical properties on the B 
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NV C 


NV D 


NV E 


NV F 


6 


35 ± 1° 


74 ± 1° 


55 ±1° 


36° ± 1° 



k° r (/xs- 1 ) 67.7 ± 3.4 63.2 ± 4.6 63.7 ± 4.5 69.1 ± 1.6 

k% (/is" 1 ) 6.4 ± 2.3 10.8 ± 4.1 9.3 ± 3.0 5.2 ± 0.8 

k° 7 (/xs- 1 ) 50.7 ± 4.4 60.7 ± 6.6 53.0 ± 5.9 48.6 ± 1.9 

k° 71 (/xs- 1 ) 0.7 ± 0.5 0.8 ± 0.6 0.9 ± 0.8 1.5 ± 0.5 

fc? 2 (/xs- 1 ) 0.6 ± 0.3 0.4 ± 0.2 0.5 ± 0.2 1.4 ± 0.2 



Table 1. Summary of the parameters extracted from the fit of time-resolved PL traces 
7Z(t, B) for four different NV defects, oriented with different angles with respect to the 
magnetic field. The indicated values are mean ± one standard deviation, calculated 
from a set of fitting parameters obtained after running the fitting procedure many 
times while varying the initial guess. 



field and its application for magnetic field imaging. Figure [4](a) shows the normalized 
PL intensity calculated using the seven-level model with the parameters of NV E as a 
function of B z = B cos9 and B± = £> sin#, varied in the range 0-150 mT. Regardless of 
the value of B z in this range, the PL steadily decreases with an increasing B±, with a 
PL drop exceeding 30% for B± > 20 mT. This effect results from a mixing of electron 
spin states both in the ground and excited levels, which enhances the probability of non- 
radiative ISC to the metastable level and reduces the PL intensity. We note that in the 
low transverse field regime, the PL intensity also exhibits the well-known sharp drops 
at B z « 51 mT and at B z w 102 mT, corresponding to spin mixing induced by a level 
anti-crossing within the excited and ground level, respectively (22,25 . In Figure [3](c), 
the normalized integrated PL signal is plotted as a function of B for three distinct NV 
defects, oriented with different angles 9 with respect to B. The data are well reproduced 
by the model (solid lines) without any free parameters. This graph highlights the fact 
that the monotonous PL decrease with B is a general feature of NV defects, occurring 
whenever B is not aligned with z to better than w 20° (see Fig. |4](a)). Such a PL drop 
can be used to discriminate between low and high transverse magnetic field regions. 
Although not fully quantitative, a PL measurement sets a lower bound to the magnetic 
field amplitude. For instance one deduces from the contour lines in Fig. [I](a) that a 
PL drop of 30% implies that the field has an amplitude B > 20 mT. This property can 
be used to develop a microscope capable of qualitatively mapping large magnetic field 
regions at the nanoscale. We note that this method could be used with any solid-state 
emitter exhibiting spin-dependent PL [34] , with a characteristic magnetic field range 
given by its zero-field splitting (Eq. ([I 

The PL drop induced by spin mixing is correlated with a reduction of the effective 
lifetime r e ff of the NV defect excited level (Fig. [3](c)), which is defined as the exponential 
decay constant that best fit a given time-resolved PL trace. Using the photophysical 
parameters of NV E, we infer the effective lifetime r e ff as a function of B z and B± by 
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Figure 4. (a) Calculated PL intensity and (b) effective lifetime r e ff as a function of 
the field components B z and Bj_, parallel and perpendicular to the NV defect axis, 
respectively. For both (a) and (b), the photophysical parameters of NV E are used 
in the seven-level model, and a few contour lines are superimposed (dotted lines), (c) 
Integrated PL intensity and (d) effective lifetime as a function of B measured for NV D, 
NV E and NV F. Error bars are 95% confidence intervals. The angles 9 corresponding 
to the three NVs considered in (c) and (d) are pictured by dashed lines in (a) and (b). 
Solid lines are the results of the calculation. 



fitting the time-resolved PL traces calculated from the seven-level model with a single 
exponential decay. As shown in Figure [4](b) , the evolution of r e ff as a function of the 
magnetic field is clearly correlated with the PL rate (Fig. |3|(a)). This feature illustrates 
that field-induced spin level mixing increases the mean probability for non radiative 
ISC transitions to the metastable level, which not only decreases the PL rate but also 
increases the decay rate from the excited level, thus leading to an overall reduction of 
the excited level lifetime. In Figure [3](d), r e ff is plotted as a function of B for three 
different defects, together with the model that introduces no free parameters. In the 
same way as for the PL rate measurement, knowing the effective lifetime sets a lower 
bound for the amplitude of the magnetic field (see contour lines in Fig. [|(b)), and can 
be used to map high off-axis magnetic field regions, as demonstrated in the next section. 




Figure 5. (a) Simplified scheme of the scanning NV defect microscope. Either a 
pulsed or a CW laser source operating at 532 nm wavelength is tightly focused at the 
end of the AFM tip through a microscope objective with a 0.9-numerical aperture and 
a 1-mm working distance (Olympus, MPLFLN100X). The PL emitted by a single 
NV defect placed at the apex of the tip is collected by the same objective and 
spectrally filtered from the remaining pump light using a dichroic mirror (DM). The 
collected light is then detected with a Hanbury Brown and Twiss (HBT) interferometer 
consisting in two single-photon counting modules (APD) placed on the output ports of 
a 50/50 beamsplitter (BS). Scanning probe magnetometry is performed by monitoring 
the magnetic-field-dependent NV defect PL while scanning the magnetic hard disk 
(HD). (b) Histogram of time delays between two consecutive single-photon detections, 
recorded in pulsed regime with the HBT detection system. The solid line is a data fit 
with multi-exponential functions. 



3. All-optical magnetic field imaging with a scanning NV defect 

Magnetic field imaging with a scanning NV defect is performed with the experimental 
setup depicted in Figure [5](a), which combines an optical confocal microscope and a 
customized tuning-fork-based atomic force microscope (Attocube Systems, CFM/AFM), 
all operating under ambient conditions. A 20-nm diamond nanocrystal hosting a single 
NV defect is first grafted at the end of the AFM tip following the procedure described in 



Refs. 19,35 . A confocal microscope placed on top of the tip allows us both to excite and 
collect the NV defect magnet ic-field-dependent PL. The unicity of the emitter placed 
at the apex of the AFM tip was checked by measuring in pulsed regime the histogram 
of time delays between two consecutive single-photon detections with a Hanbury Brown 
and Twiss interferometer. As shown in Figure [5](b), a strong reduction of coincidences 
is observed around zero delay, which is the signature that a single NV defect is attached 
to the tip [36] . 

As a test sample, we used a piece of a commercial magnetic hard disk composed 
of 450 nm x 55 nm bits with random in-plane magnetization (Fig. [5](a)). Two adjacent 
magnetic bits with opposite magnetizations produce a stray field coming in or out at the 
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Figure 6. (a) PL map recorded by scanning the magnetic hard disk while operating 
the AFM in tapping mode, with CW optical excitation. The image corresponds to 
300 x 300 pixels, with a 8-nm pixel size and a 20-ms acquisition time per pixel. The 
white dashed lines indicate the tracks of the magnetic hard disk, (b) PL decay traces 
recorded with the AFM tip standing either on a bright or on a dark area. Solid lines are 
data fits using single exponential decays. (c,d) PL (c) and effective lifetime (d) maps 
obtained under pulsed laser excitation. The scanned area (39 x 30 pixels) corresponds 
to the dashed red square depicted in (a) with a 3-s acquisition time per pixel. 

boundary between the two bits. Imaging such a stray field is performed by recording 
the NV defect PL intensity while scanning the magnetic hard disk. In this experiment, 
the AFM is operated in tapping mode in order to maintain the mean probe-to-sample 
distance to a constant value. As shown in Figure [6](a), the PL image reveals an array of 
dark areas that correspond to regions of high transverse magnetic field. Such dark areas 
indicate magnetization reversals, very much like in magnetic force microscopy (MFM) 
measurements with the difference that in addition MFM distinguishes between fields 
coming in or out of the sample, corresponding to repulsive or attractive forces. The 
distance between two consecutive dark lines is measured to be larger than 100 nm in 
Figure [6](a). For instance, it is approximately 110 nm in the zoom shown in Figure [6](c), 
which corresponds to the size of two magnetic bits. We therefore believe that in this 
particular experiment, at least two consecutive bits with the same magnetization are 
needed in order to build up a stray field that is strong enough to induce a detectable 
PL quenching. 



Magnetic-field- dependent photodynamics of single NV defects in diamond: Application to qualitative all-opti 



As discussed in Section [2j the field-induced PL quenching is associated with 
an overall reduction of the NV defect excited-level lifetime. This is illustrated in 
Figure [6](b) , which shows PL decay traces measured when the tip stands either on a 
bright area (low magnetic field) or on a dark area (high transverse magnetic field). Single 
exponential data fits give effective lifetimes Tbright — 18.2 ± 0.3 ns and Td ar k = 15.1 ± 0.2 
ns, corresponding to a reduction of 17%. We note that the lifetime of the NV defect is 
larger in nanocrystals compared to bulk measurements reported in Figure [|(b)), owing 
to change of the surrounding refractive index [37] . By recording PL decay traces while 
scanning the sample, magnetic field mapping can also be achieved through fluorescence 
lifetime imaging (FLIM), as shown in Figure [6](d)). As expected, FLIM and PL images 
are clearly correlated: the regions of low effective lifetime closely match those of low PL 
signal, in correspondence to regions of high off- axis magnetic field (Fig. [6](c)-(d)). 

The maximum PL contrast in Fig. [6](b) is about 30%, which is close to what was 
measured for NV defects in bulk diamond. As discussed in Section [2j measuring the 
normalized PL intensity can in principle give access to a lower bound for B provided that 
one knows all relevant parameters of the experiment, including the transition rates of the 
NV defect employed. However, those are affected by the electromagnetic environment 
of the NV defect - proximity to the AFM tip, to the metallic sample, etc. - and are 
therefore not readily accessible. Assuming that the general behavior of Figure [I](a) is 
however still approximately valid here, one can state that the dark regions in Figure |6](a) 
correspond to a magnetic field exceeding 10 mT. This would be consistent with the stray 
field expected for such a magnetic hard disk at a distance of a few tens of nanometers 1 19 1 . 
Although not quantitative, all optical magnetic field imaging with a scanning NV defect 
is thus a relatively simple way to map regions of magnetic field larger than a few tens of 
mT at the nanoscale. Furthermore, unlike MFM whose magnetic tip is likely to perturb 
the magnetization of the studied sample, the NV defect probe introduces no magnetic 
back action to the sample at first order. 

Finally, we note that the ability to perform fluorescence lifetime imaging (FLIM) 
with a scanning NV defect opens up a large range of possibilities in the field of 
nanophotonics. Indeed, mapping the local density of electromagnetic states (LDOS) 
above non-magnetic samples is in reach, with the potential for vectorial LDOS mapping 
if one uses NV defects with different orientations. Compared to recent demonstrations 
of scanning FLIM (20|[2l], advantages of NV defects are their perfect photostability as 
well as their atomic-sized spatial resolution. In this context, this work is a starting-point 
for understanding and quantifying the relation between the NV defect dynamics and the 
LDOS, a necessary step for future LDOS mapping experiments. 

4. Conclusion 

In this article the optical properties of NV defects in diamond have been studied as a 
function of the external magnetic field. We performed time-resolved measurements and 
developed a seven-level model that accounts for the decreased ESR contrast, PL intensity 
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and effective lifetime when a transverse magnetic field is applied. We demonstrated an 
application of those effects to nanoscale magnetic imaging. Using a scanning NV defect 
microscope, we mapped the stray field of a magnetic hard disk by recording either PL 
or lifetime images. This all-optical method for high magnetic field imaging might be 
of interest in the field of nanomagnetism, where samples producing fields in excess of 
several tens of milliteslas are typical. 
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